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Abstract 

We present a Mathematica notebook allowing for the symbolic calculation of the 3x3 
dielectric tensor of a electron-beam plasma system in the fluid approximation. Cal- 
culation is detailed for a cold relativistic electron beam entering a cold magnetized 
plasma, and for arbitrarily oriented wave vectors. We show how one can elaborate 
on this example to account for temperatures, arbitrarily oriented magnetic field or 
a different kind of plasma. 
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1 Introduction 



The calculation of the dielectric tensor of a beam plasma system is a recurrent 
problem in plasma physics. Many efforts have been dedicated recently to such 
issue because of the Fast Ignition Scenario for inertial thermonuclear fusion 
pl2] . According to this view, the Deuterium Tritium target is first compressed 
by means of some driver. Then, the compressed fuel is ignited by a relativistic 
electron beam generated by a petawatt laser shot. Such scenario implies there- 
fore the interaction of a relativistic electron beam with a plasma. This kind of 
interaction, and its magnetized counterpart, is also relevant to astrophysics, 
in particular when investigating the relativistic jets of microquasars [3], active 
galactic nuclei [I], gamma ray burst production scenarios [5 J or pulsar winds 
[6]. Theoretical works on these subjects are usually focused on the instabilities 
of the system. Although many of them demands a kinetic treatment to be fully 
described, the fluid equations can set some very relevant guidelines, especially 
when the system is not too hot. Furthermore, it has been known for long 
that in the relativistic regime, instabilities with arbitrarily orientated wave 
vectors may be essential [Tll8ll9lfT0][TT] . One can therefore figure out how some 
refined kinetic theory may lead to almost unsolvable calculations whereas the 
fluid formalism is still tractable. For example, a detailed description of the 
collisional filamentation instability (k _L beam) including the movement of 
the background ions plasma, and accounting for temperatures, was first per- 
formed through the fluid equations [12]. The very same equations were used 
to explore the growth rate of unstable modes with arbitrarily oriented wave 
vectors (with respect to the beam) when a relativistic electron beam enters 
a plasma [71l9lll0llllj . The results were found crucial as it was demonstrated 
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that the fastest growing modes were indeed found for obliquely propagating 
waves, and the kinetic counterpart of these models has only been considered 
very recently [T51ll4|ll5j . As far as the magnetized case is concerned, the kinetic 
formalism has been thoroughly investigated for wave vectors parallel and nor- 
mal to the beam [T6lll7|ll8] . But the unstable oblique modes, which once again 
turn to be the most unstable in many cases, could only be explored through 
the fluid formalism [8]. 

It has been demonstrated that the fluid equations yield the same first order 
temperature corrections than the kinetic theory for oblique modes, and the 
roles of both beam and plasma parallel and perpendicular temperatures are 
retrieved [TJJ]. The fluid approximation is thus definitely a tool of paramount 
importance to deal with beam plasma instabilities. Additionally, it generally 
yields a polynomial dispersion equation for which numerical resolution is im- 
mediate. Nevertheless, even the fluid tensor can be analytically involved when 
considering arbitrarily oriented wave vectors, a guiding magnetic field, tem- 
peratures, and so on [20]. Indeed, on can think about any model based on 
whether the system is relativistic or not, collisional or not, magnetized or not, 
hot or cold. . . Most of these models have not been implemented yet, and each 
one should leave a quite complicated dielectric tensor. 

This is why a Mathematica notebook has been developed which allows for the 
symbolic calculation of the fluid tensor, once the parameters of the system 
have been set. The basic system we study here is a cold relativistic electron 
beam entering a cold magnetized plasma with return current. As the reader 
shall check, the notebook is very easy to adapt the different scenarios (ion 
beam, temperatures, pair plasma...). The paper is structured as follow: we 
start introducing the theory leading to the fluid dielectric tensor in section [2j 
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The Mathematica notebook is then explained step by step in section [3j and 
we show how it can be modified to include temperatures or collisions before 
the comments and conclusion section. 



2 Theory 

We consider a beam of density n b , velocity V 6 and relativistic factor 7 b = 
1/(1 — V b 2 /c 2 ) entering a plasma of density n p . Ions from the plasma are con- 
sidered as a fixed neutralizing background, and an electron plasma return 
current flows at velocity V p such as n p ~V p = n^V;,. The system is thus charge 
and current neutralized. We do not make any assumptions on the ratio nb/n p 
so that the return current can turn relativistic for beam densities approaching, 
or even equalling, the plasma one. We set the z axis along the beam velocity 
and align the static magnetic field along this very axis. The wave vector inves- 
tigated lies in the (x, z) plan without loss of generality [Sj, and we define the 
angle 9 between k and V& || Bo || z through k z = kcosO and k x = ksin6. The 
dielectric tensor of the system is obtained starting with the fluid equations for 
each species j = p for plasma electrons and j = b for the beam ones, 



where pj = 7jmvj, m the electron mass and q < its charge. The equations 
are then linearized according to a standard procedure [8], assuming small 
variations of the variables according to exp(ik • r — iut). With the subscripts 
and 1 denoting the equilibrium and perturbed quantities respectively, the 
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linearized conservation equation ([I]) yields 

%i = n j0 , (3) 

u — k-Vjo 



and the force equation (j2J) gives, 

im 7j (k • v i0 - w) U jx + ^-(v i0 • Vji)vjoJ 
= ^E 1+ (Vj0 + Vjl)X c Bo + V ^ XBl ), (4) 

where = —1. Through Maxwell- Faraday equations, the field Bi is then 
replaced by (c/u)k x Ei so that the perturbed velocities v.,-1 can be explained 
in terms of E x resolving the tensorial equations (jl]). Once the velocities are 
obtained, the perturbed densities can be expressed in terms of the electric field 
using Eqs. ([3]). Finally, the linear expression of the current is found in terms 
of Ei through, 



J = q n jQ Vji + UjtVjo, (5) 

j=P,b 



and the system is closed combining Maxwell Faraday and Maxwell Ampere 
equations, 

r 2 4-77T 

— kx(kxE 1 )+E 1 + J = 0. (6) 

UJ Z CO 



Inserting the current expression from Eq. (jSJ) into Eq. ([6]) yields an equation 
of the kind T(Ei) = 0, and the dispersion equation reads detT = 0. 

The Mathematica notebook we describe in the next section performs a sym- 
bolic computation of the tensor T and the dispersion equation detT = 0, in 
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terms of the usual [21] reduced variables of the problem 

Z = , x = — , a = — , p=—, Ub = — , (7) 

UJp UJp Tip C LOp 



where <Ji = 47m p q 2 /m is the electron plasma frequency and u b = \q\B$/mc 
the electron cyclotron frequency. 



3 Mathematica implementation 

For the most part, Mathematica is used to solve the tensorial equations (jlj) for 
Vji and extract the tensors T from Eqs. (|5|6|) . We start declaring the variables 
corresponding to the wave vector, the electric field, the beam and plasma drift 
velocities and the magnetic field, 

In[l]:= k = {kx, 0, kz}; El = {Elx, Ely, Elz}; VOb = {0, 0, Vb}; 
VOp = {0, 0, Vp}; B0={0, 0, m c cob/q}; Bl = c Cross[k, 
vbl = {vblx, vbly, vblz}; vpl = {vplx, vply, vplz}; 

Note that Maxwell Faraday's equation is already implemented in the definition 
of Bl. The wave vector has no component along the y axis and the beam and 
plasma drift velocities only have one along the z axis. The guiding magnetic 
field is set along z and defined in terms of the cyclotron frequency u b . This 
will be useful later when introducing the dimensionless parameters (j7]). 

We then have Mathematica solve Eqs. (jlj) for the beam and the plasma. The 
left hand side of the equation is not as simple as in the non-relativistic case 
because the 7 factors of the beam and the plasma modify the linearization 
procedure. We write this part of the equations in a tensorial form in Mathe- 
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matica defining the tensors Mp and Mb such as "left hand side"=Mj 1 .vjl 
with, 

In[2]:= Mb={{ 7b( J fczyfe) ,0,0},{0, 7 ^_ i fc ^ b) ,0},{0,0, 7b3( J fczyfe) }}; 

In[3]:= Mp={{ 7p( J fczyp) ,0,0},{0, 7p( J fc ^ p) ,0},{0,0, 7pa( J^ p) }}; 

where i 2 = —1. The reader will notice that relativistic effects are more pro- 
nounced in the beam direction due to the 7 3 factors in the zz component. We 
now have Mathematica solve the tensorial Eqs. (j3J). For better clarity, we first 
define them 

Jn#/:=EqVb=vbl-Dot[Mb,£ (El + Cross[VOb+vbl,BO] ) ]_ Dot [Mb,± ( Cross [VOb.BlJ j]. 

/n/5/. =EqVp=vpl-Dot[Mp,^ (El + ^[VOp+vpi,Bo] ) ]-Dot [Mp, ^ ( CroSB ^ 0p ' B1] )]; 
before we solve them, 

In[6j .=Vbl=FullSimplify[vbl/. Solve[EqVb==0,vbl][[l]]]; 

Jn/7/:=Vpl=FullSimplify[vpl/. Solve[EqVp==0,vpl][[l]]]; 

Note that the Vb's, with capital "V", store the solutions of the equations 
whereas the vb's are the variables. This is why the Vb's do not need to be 
defined at the beginning (see In[lJ) of the notebook; they are implicitly defined 
here. 

Now that we have the values of the perturbed velocities, we can derive the 
perturbed densities from Eqs. d3J), 

/n/g/.=Nbl=FullSimplify[^ 

/n^.=Npl=FullSimplif y [^pp^^gg^]; 
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Here again, we prepare the introduction of the reduced variables (J7J) by ex- 
pressing the equilibrium beam and plasma electronic densities in terms of the 
beam and plasma electronic frequencies. 

We can now have Mathematica calculate the current according to Eq. fl5j), 
In[10]:= 

J=FullSimplify[q(c,;pp 2 i ^Vpl + cjpb 2 ^ Vbl+NplVOp + NblVOb)]; 

We now have the symbolic expression of the current J. In order to find the 
tensor T yielding the dispersion equation, we need to explain first the current 
tensor. This is performed through, 

In[ll]:=M= 

( \ 
Coefficient [ J [[1]], Elx] Coefficient [ J [[1]], Ely] Coefficient [J [[1]], Elz] 

Coefficient [J [[2]], Elx] Coefficient [J [[2]], Ely] Coefficient [J [[2]], Elz] 5 

Coefficient [J [[3]], Elx] Coefficient [J [[3]], Ely] Coefficient[J[[3]], Elz] 



which just extract the tensor elements from the expression of J. We now turn 
to Eq. (jBJ) where we explain the tensor elements of the quantity c 2 k x (k x 
E x ) +w 2 Ei, 

In[12]:=M0=c 2 Cross[k,Cross[k,El]]+cj 2 El ; 
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In[13]:=Ml = 

( \ 
Coefficient[MO[[l]], Elx] Coefficient[MO[[l]], Ely] Coefficient[MO[[l]], Elz] 

Coefflcient[M0[[2]],Elx] Coefficient[M0[[2]], Ely] Co efficient [MO [[2]], Elz] 5 

^ Coefficient[M0[[3]], Elx] Coefficient[M0[[3]], Ely] Co efficient [MO [[3]], Elz] 

We can finally express the tensor T defined by T(E)=0 as 
In[14j :=T=Ml+4 i n uj M; 

At this stage of the notebook, we could take the determinant of the tensor 
to obtain the dispersion equation. Let us first introduce the dimensionless 
variables ([7]) through, 

In[l 5] :=T=T /. {Vp -> -a Vb, kz -> upp Zz/Vb, kx -> cupp Zx/Vb, 
c<jpb 2 — > a cupp 2 , c<j — > x c<jpp, uh — > fib tupp}; 



and, 



T=T /. {Vb -> /3 c} 



Mathematica leaves here some cupp's 



which should simplify between each oth- 



ers. It is enough to perform 



In[17]:=T=T /. {upp - 1}; 



and a simple 



/n/lS/.-=MatrixForm[FullSimplify[T]] 



displays the result. The dispersion equation of the system is eventually ob- 
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tained through 
/n/15/.-=DisperEq=Det[T] 

The notebook evaluation takes 1 minute on a Laptop running a 1.5 GHz 
Pentium Centrino under Windows XP Pro. This delay can be shortened down 
to 10 seconds by suppressing all the FullSimplify routines while leaving a 
Simplify[T] in entry 18, but the final result is much less concise and readable. 



4 Comments and Conclusion 

In this paper, we have described a Mathematica notebook performing the 
symbolic evaluation of the dielectric tensor of a beam plasma system. Start- 
ing from the linearized fluid equations, the notebook expresses the dielectric 
tensor, and eventually the dispersion equation, is terms of some usual dimen- 
sionless parameters. This notebook has been so far applied to the treatment 
of the temperature dependant non magnetized and magnetized problems (see 
Refs [T9"f2"0"] ). Indeed, the procedure is very easy to adapt to different settings. 

When including beam or plasma temperatures, one adds a pressure term 
—VPj/rij on the right hand side of the force equations (J2J). Setting then 
VPi = ?>kBT,jVni [12.22] if dealing only with electron motion, one only needs 
to add to the notebook entries 4 an d 5 the terms (i 2 =-l) 

o: rp- i Dot[k,vjl] 
"«« X J * w _Dot[k,VOj]' 

where j=p for the plasma, and b for the beam. When considering anisotropic 
temperatures [19], one just needs to define a temperature tensor Tj for each 
species j, and replace the scalar product Tj k by the tensorial one Dot[Tj,k] 
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in both entries. Of course, a correct treatment of electromagnetic instabilities 
generally requires a kinetic formalism instead of a fluid one. However, kinetic 
calculations cannot be systematically entrusted to Mathematica, as is the case 
here. The reason why is that the relativistic factors 7 encountered in the 
kinetic quadratures are coupling the integrations along the three components 
of the momentum. According to the distribution functions considered, the 
quadratures may be calculable through some ad hoc change of variables, if 
they can be calculated at all. At any rate, the process cannot be systematized 
enough for Mathematica to handle it. 



As far as the magnetic field is concerned, its direction can be changed from 
entry 1 without any modification of the next entries. When dealing with the 
motion of ions, or even with one of these pair plasmas involved in the pulsar 
problems [23], one just need to modify the conservation and force equations 
according to the properties of the species investigated. It is even possible to 
add more equations to account for more species because the resolution involves 
only the force and the conservation equations of one single specie at a time 
before the perturbed quantities merge together in entry 10 to compute the 
current J. 



The notebook can thus be easily adapted to different settings and allows for a 
quick symbolic calculation of the dielectric tensor and the dispersion equation, 
even for an elaborated fluid model. 
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